Consensus in networks of mobile communicating agents 
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Populations of mobile and communicating agents describe a vast array of technological and natural 
systems, ranging from sensor networks to animal groups. Here, we investigate how a group-level 
agreement may emerge in the continuously evolving network dened by the local interactions of 
the moving individuals. We adopt a general scheme of motion in two dimensions and we let the 
individuals interact through the minimal naming game, a prototypical scheme to investigate social 
consensus. We distinguish different regimes of convergence determined by the emission range of the 
agents and by their mobility, and we identify the corresponding scaling behaviors of the consensus 
time. In the same way, we rationalize also the behavior of the maximum memory used during 
the convergence process, which determines the minimum cognitive/storage capacity needed by the 
individuals. Overall, we believe that the simple and general model presented in this paper can 
represent a helpful reference for a better understanding of the behavior of populations of mobile 
agents. 
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I. INTRODUCTION 

Autonomous mobile and communicating agents pro- 
vide extremely efficient solutions to a wide range of tech- 
nological problems by guaranteeing robustness, flexibil- 
ity, and dynamic adaptability A typical case is that 
of a population of robots that have to explore an un- 
known environment, and cope with situations that by 
definition can not be foreseen in advance 0, Q . For ex- 
ample, robots could have to negotiate a common lexicon 
to name different places of the environment they are sur- 
veying, and then use this shared linguistic knowledge to 
carry out goal-directed behavior 0, IB] ■ In the same way, 
also the performances of sensor networks can be en- 
hanced by the introduction of mobile agents 0, Q ■ In 
natural systems, on the other hand, mobile populations 
that coordinate through chemical or audible signals are 
obviously widespread, ranging from cell populations Q 
to animal groups [lOj . 

In all of these cases, mobil e agent s locally broadcast 
their signal to nearby nodes [lll - ll3| . Thus, communi- 
cation takes place on a continuously evolving network 
whose properties are determined by such parameters as 
the emission range of the individuals, their mobility, or 
their density. Network theory (l3 - [T8l | is therefore the nat- 
ural framework to investigate the emerging population- 
scale properties of the system. However, previous re- 
search has so far focused mainly on static random ge- 
ometric networks [ill, [Hi, or on the opposite case of 
rapidly changing structures [l^, [13] ■ Only very recently 
the more general case of time-dependent networks has 
been fully addressed, for the specific case of the synchro- 
nization of mobile oscillators |21| . 

In this paper we address the fundamental problem of 
social consensus. To this purpose, we model individuals 



that move in a 2-dimensional plane and have to agree 
on a given convention by performing standard language 
games [12, [l^- For example, they might be in need of 
creating or selecting autonomously a key for encrypted 
communication |l2|, or to independently elect a leader 



[2J| . We then study how different parameters of the mo- 
bile agents, in particular their velocity and their commu- 
nication range, affect the overall agreement process. We 
are able to identify different regimes ruling the consensus 
dynamics, and we rationalize our findings by considering 
the properties of the communication networks. 



II. THE MODEL 

We model the mobility of the agents according to the 
general scheme put forth in [ll|. A population of N 
individuals move in a 2-dimensional L x L box with pe- 
riodic boundary conditions. The velocity of each agent 
is V. The angle of the ith agent's motion is £,i{tk) € 
[0, 27r], and it changes randomly at discrete time steps tk 
(tk+i — tk = tm). Thus, the evolution in time of the ith 
agent's position is: 



Xt{tk + At) = Xt{tk) + V cos ^i{tk) At mod L (1) 
yi{tk + At) = y,{tk)+v sin ^i{tk) At modi, (2) 

where A< < rj\/. The motion of the individuals is there- 
fore diffusive, the diffusion coefficient being D ^ v^tm- 

The agents play the minimal naming game (NG) 
[25! . [26l | implemented with local broadcasting without 
feedback [l3[. Each agent is characterized by an inven- 
tory of words (or "conventions" , "opinions" , "forms" or 
"states"). At the beginning all inventories are empty. At 
discrete time steps of duration an agent is randomly 
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selected as speaker. She selects randomly a word from 
her inventory, and transmits it to all the agents within 
a distance d from her position (if the inventory is empty 
she invents a brand new word and stores it into her own 
inventory, before broadcasting it). Each receiving agent 
updates her state depending on whether her inventory 
contains or does not contain the transmitted word. In 
the first case, the agent deletes all the competing syn- 
onyms and keep only that word into her inventory. In 
the latter case, on the other hand, she adds the new re- 
ceived word to her inventory. The speaker receives no 
feedback about her emission, and consequently does not 
modify her inventory Figure [1] summarizes the rules of 
the model. 

For simplicity, in this paper we choose At = tm = 
T5 = 1. This means that in a time step (i) all agents move 
in straight line, (ii) all agents are reassigned a random 
angle, and (iii) one agent broadcasts to her neighbors. 



III. PATHS TO CONSENSUS 

The NG is an ordering process. An initially disordered 
configuration ends up in a consensus, ordered, state in 
which everybody has the same unique word [1^ H^] . The 
consensus (or "convergence") time, tconv, is therefore a 
crucial quantity. Also important is the maximum num- 
ber of words agents have to store, M, which describes 
the global amount of memory needed by the system to 
reach a consensus. Previous studies on quenched graphs 
have shown that both quantities depend dramatically on 
the topology of the social network describing the possible 
interactions between the individuals [28| . Consequently, 
to investigate the properties of a mobile population it is 
convenient to focus on the properties of the static net- 
work describing the instantaneous communications of the 
agents. This is the graph that is obtained, at any time, 
by drawing an undirected link between any two agents 
that are closer than the emission range d. Recalling that 
the average degree (i.e., the average number of neighbors 
of a randomly selected node) of the network is simply 
(k) = irNd'^/L'^, these values of d identify different sce- 
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is the range above which the average 



degree is larger than 1, so that every emission is 
received on average by some agent. 



dr = d 



1 is the critical value for a percolation 
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transition, yielding a giant component of size N 

dj^/2 = c^{fc)=Af/2 is the point where every communi- 
cation involves on average the majority of the pop- 
ulation. 



• dr, 



= d 



{k)=N 
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FIG. 1: The model. N agents (circles) move with velocity 
V and randomly assigned angles in a box of size L, with 
periodic boundary conditions. At each time step one of the 
agents is chosen as a speaker (black circle) and emits a word 
randomly extracted from her inventory [A,'m figure). All the 
agents within a distance d (gray circles) receive the word, and 
update their inventories as shown in the schematic represen- 
tation below the box. If an agent already knows that word, 
she deletes all the competing synonyms in her inventory, oth- 
erwise she simply adds the new word to it. No feedback is 
provided to the speaker, whose inventory is not altered. 



Of course, it holds di < dc < d]\[/2 < dmax- In this paper 
we fix L = 200 and N = 100, unless where explicitly 
stated, so that di ~ II. 3, dc ~ 24.1, djv/2 — 79.8 and 
dmax — 141.4. Figure [5] shows the dependence of the 
average degree and the size of the giant component on 
the parameter d for this choice of L and N. 

For a qualitative partitioning of the observed phe- 
nomenology in terms of distinct regimes, it is convenient 
to consider two timescales. One describes the stability of 
a cluster of agents, and the other accounts for the time 
over which a consensus is reached within the same cluster 
[3^. Their ratio 77 assesses therefore the impact of local, 
intra-clustcr, activity on global, inter-cluster, dynamics. 
For the robustness of a cluster, we consider the average 
number of timesteps needed by an individual to leave a 
group of size n(d) (an increasing function of d for d < dc), 



that scales as ti ~ n{d)/v'^ [2l|. For the within-cluster 



nected network. It holds dr, 



/LV2. 



average consensus time, on the other hand, we note that 
it can be treated as independent from the cluster size 



FIG. 2: (Color online) Properties of the static network. At 
d > di every emission is lieard on average by at least an in- 
dividual of the population, while at d = djv/2 a majority of 
agents listens to each communication act. At d = dc the giant 
connected component ("gcc") is formed by N nodes (in real- 
ity, due to finite size effects, this happens for a slightly larger 
emission range). At dmax the network is fully connected. 

n{d), tconv ^ const., both when the considered groups 
are densely connected and when they are very small. The 
reason is that in the first case the broadcasting rule brings 
about a very fast consensus time, which becomes instan- 
taneous in fully connected graphs [l^ , while in the latter 
case consensus is quick simply because just a few agents 
have to agree [l3l l25l|. This approximation is appropri- 
ate for our purposes since we aim at defining a qualita- 
tive index able to discriminate between extreme regimes. 
Moreover, as we shall see, it is further validated by the re- 
sults discussed in the following of the paper (Sec. IIII Cp . 
Therefore, for the ratio r] between the two timescales it 
holds: 




FIG. 3: (Color online) Top: Consensus time as a function of d 
and V (all axes report the logjQ of the respective quantities). 
Bottom: t^onv as a function of d for different for different 
values of the agents velocity v. Dotted vertical lines represent 
di, dc and djv/2. 



which is obviously an increasing function of d as d < dc. 

A. Consensus time 

Figure[3]shows the behavior of the consensus time tconv 
as a function of the emission range d, and for different 
values of the agents' velocity v. The consensus is fast for 
large values of d (and becomes instantaneous as soon as 
d = dmax , when everybody receives the word transmitted 
by the first speaker), but it increases for shorter ranges, 
in a way that crucially depends upon the v parameter. 
We can identify three regimes. 

(1) 77 <C 1 holds for small d and large w, and implies 
a rapidly evolving network. Agents continuously change 
their neighbors, and hence the partners of their commu- 
nication acts, pretty much as in the case of an annealed 



network. Thus, consensus emerges through global agree- 
ment at the system size level, after the agents have cor- 
related their inventories so as to allow for successful com- 
munication to take place (soj . As d < di the behavior 
tconv ~ = is observed (Figure SI a), describ- 

ing the existence of empty communication acts (unheard 
emissions) when on average each node has less than one 
neighbor. 

(2) 77 > 1 and d < dc, on the other hand, implies 
smaller velocities. In this case, small and isolated clus- 
ters of agents locally reach an agreement on different 
conventions. Global consensus emerges at a later time 
through the competition between these words, in a sit- 
uation reminiscent of what happens in low-dimensional 
lattices (Slj . The intra-cluster movements determine 
the leading timescale, implying a scaling of the form 
tconv ^l/v^ (Figure H b and c). 

(3) 77 > 1 and d ^ dc, finally, identify a scenario in 
which the whole population forms a single connected clus- 
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FIG. 4: (Color online) Rescaling of the consensus time, d < 
dc- When 77 <C 1, the d^tconv ~ const, and curves for large 
velocities collapse (a). For > 1, on the other hand, curves 
for different, and small values, of v collapse as the consensus 
time is rescaled as v^tconv ~ const (b). This behavior is 
observed also for values of d close to dc provided that small 
enough velocities are considered (c). In all panels horizontal 
dashed lines represent a constant behavior, and serve as a 
guide for the eye. 

ter, describing a random geometric graph. In [l^ Lu 
and coworkers showed numerically that tconv ^ l/{k)^'^ 
( "when (fc) ^ N") for static random geometric graphs 
[341 . Accordingly, in Figure [5] we observe the behavior 
tconv ^ l/d^'^, which as expected degrades before dmax, 
where (fc) = N — 1. It may be further noted that, as 
d > (ijv/2, the curves for different v behave identically 
(Figure [5l inset). This is due to the fact that here the 
first speaker transmits her word to an absolute majority 
of the agents, which on their turn drive the system to 
consensus very rapidly thanks to the fact that the NG is 
a drift-driven process [2^ . 



B. Memory usage 

The agents get to know different words at the same 
time during the process that eventually leads them to a 
consensus, and the NG rules do not fix a limit to the 
size to their inventory. Therefore, it is important to look 
at the maximum memory consumption that the popula- 
tion experiences during the whole process, correspond- 
ing to the maximum number of words, M, present in the 
system, i.e., to the sum of the inventory sizes of the N 
agents. Gaining quantitative insights into the behavior 
of this quantity is more difficult than for tconv , but some 
hints can be gained from numerical investigations. Fig- 
ure [5] shows the average maximum memory per agent as 
a function of d, and for different values of v. 

Again, to understand what goes on in the case of mo- 
bile agents it is helpful to recall the results obtained in 
static networks. It turns out that a finite connectivity im- 



FIG. 5: (Color online) Rescaling of the consensus time, 
d > dc- For large 77 > 1 and large d curves for different values 
of 1; behave as tconv ~ as far as d (and hence (k)) is not 

too large. The behavior is better observed for smaller veloci- 
ties, since large values of v reduce the value of rj. As d > df^/2 
the different curves collapse, and v becomes an irrelevant pa- 
rameter, as shown also in the inset with non-rescaled abscis- 
sas. Horizontal dashed lines represent a constant behavior, 
and serve as a guide for the eye. 



plies finite memory requirements [28[, while a fully con- 
nected graph would require infinite inventories (in the 
thermodynamic limit) [25|. In general, a larger average 
degree requires a bigger memory effort for the agents |28| . 
Of course, however, the broadcasting rule implies an im- 
mediate consensus on fully connected graphs, entailing 
a minimal amount of memory. Moreover, ordered low- 
dimensional lattices determine an extremely reduced use 
of memory, since convergence is reached through the com- 
petition of different clusters of agents who have reached a 
local consensus and therefore store one word only j30ll3ll |. 

In the light of these results, it is possible to rationalize 
the findings presented in Figure [51 As expected, when 
the dominant process is global agreement (77 <C 1, small d 
and large v), more memory is needed [region (1), above]. 
Every agent is exposed to a large number of words due 
to the high mobility of the population. As the velocity 
is decreased, and the system enters the phase of isolated 
clusters, on the other hand, consensus is reached with a 
smaller memory demand, due to the early onset of regions 
of local agreement [region (2)]. As the emission range is 
increased, and one single connected cluster emerges, fi- 
nally, curves for different values of v become more and 
more similar, and collapse as d > diy^2 (Figure [51 inset), 
as observed for tconv [region (3)]. Curves for small v ex- 
hibit a peak somewhere in the region is dc < d < djv/2 
(more precisely for 50 < d < 60 in Figure [6]), separating 
the memory-favorable cases of a fully connected graph 
and a collection of isolated clusters. However, as d is 
further decreased, the memory requirements grow again. 
Qualitatively, this behavior can be understood consid- 
ering that when d < di, local clusters are in general 
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FIG. 6: (Color online) Memory usage. The average maximum 
memory used by each agent is plotted for different values of 
the velocity v. While for large mobility the individual storage 
capacity monotonously increases as d is reduced, small veloc- 
ity induces a more complex behavior. In this case, the curves 
exhibit a minimum for d ~ di and a maximum for d ~ 60, for 
the usual choice of the parameters. For d > djv/2 curves for 
different velocities collapse (see Inset). 
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FIG. 7: (Color online) Role of the population size (at fixed L). 
Top: The consensus time reaches a plateau at large in the 
case of small velocities (left panel), while it grows very weakly 
for larger mobilities (right panel). Bottom; The maximum 
memory per agent required during the process is constant at 
large TV for both high and low mobility rates. 



tree-like, so that more words can coexist in each of them 
before a local agreement emerges. 



C. Role of the population density 

In finite static networks the global consensus is reached 
only if d > dc, when all the nodes belong to the same 



unique giant connected component. Previous studies 
[l^ . conducted adopting a sHghtly modified NG pro- 
tocol, showed a A^", with a ~ 2.1, dependence of the 
consensus time on the population size N in random geo- 
metric graphs with fixed average connectivity (and con- 
sequently varying L), and, as mentioned above, an im- 
portant role also for the average degree, tconv ^ (fc)^^-^ 
(as far as fc ^ N) for fixed L. The scaling of tconv with 
the population size agrees well with a theoretical argu- 
ment put forth in [l2j , which is in its turn very similar to 
the analysis of the NG in low-dimensional lattices [3l| . 
The crucial point is that different clusters of local con- 
sensus emerge rapidly, and the global agreement is the 
result of a cluster-cluster competition. The presence of 
mobility makes this whole argument break down, and 
the increased population mixing yields a faster scaling of 
the consensus (data not shown), i n ag reement to what is 
observed in small- world networks |l2l l32j. 

Crucially, on the other hand, mobility guarantees that 
the consensus is always reached, at least asymptotically, 
independently from the emission range of the individuals. 
However, it is important noting that population size and 
average degree are intimately connected in the framework 
we are addressing, where the box side L is kept constant. 
Indeed, as discussed above, it holds (fc) cx N, and there- 
fore 



N^oo 



0, 



(4) 



and (ijv/2 being the only characteristic length surviving 
in the thermodynamic limit. The "small" range regimes 
we have described would therefore vanish for very large 
population size. 

Concretely, the case of very large population densities 
is quite unrealistic from the point of view of any applica- 
tion involving mobile individuals. However, for the sake 
of completeness, we report in Figure [7] the results of nu- 
merical investigations for different values of d and v. For 
small velocities (left column), very small population sizes 
are not efficient since the individuals waste time in find- 
ing each other. As d grows, however, tconv saturates. 
After a certain threshold, the broadcasting rule makes 
the actual number of agents present in the population 
irrelevant. A similar behavior is observed for the case of 
a larger mobility of the individuals (right column). Here, 
however, smaller population sizes converge slightly faster 
since agents get more easily in contact, and the dynamic 
of the NG favors a smaller number of competing conven- 
tions. In addition, the augmented population mixing is 
responsible for the fact that, (i) the consensus time docs 
not reach a plateau but keeps increasing very slowly with 
the system size, as tconv ^ N", with a < 0.05 and (ii) 
the memory consumption is larger, in agreement with 
what is observed in the case of static networks for the 
NG with pairwise interactions [l^. Finally, it is worth 
noting that the very weak dependence of the consensus 
time on N validates, a posteriori^ the assumption on the 
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constant behavior of tconv that yielded to Eq. [3] 

IV. CONCLUSION 

In this paper we have analyzed the consensus problem 
in a population of autonomous mobile agents. We have 
focused on the crucial case of a self-organized agreement 
process, to be established without any central control nor 
coordination. Agents move with the same velocity and 
different, randomly changing, angles in a two dimensional 
space, and communicate through the NG protocol, lo- 
cally broadcasting in an circle of radius d. We have shown 
that different characteristic emission ranges exist, defin- 
ing, together with the mobility rates of the agents, the 
boundaries between different consensus regimes. In par- 
ticular, we have highlighted three main mechanisms rul- 
ing the onset of consensus. If the emission range is small, 
a rapidly mixing population will undergo a global agree- 
ment process, while a slower mobility will bring about the 
appearance of isolated clusters in which a local consensus 
on different conventions forms rapidly, the final consensus 
resulting from the competitions of these clusters. Finally, 
large emission ranges establish a single connected cluster 
where the static limit of random geometric networks is 
recovered. Accordingly, we have pointed out the scaling 
relations of the consensus time in each region and ratio- 
nalized the memory needs of the agents. Finally we have 
considered the role of the population density, showing 
that, due to the broadcasting rule, even unrealistically 
high densities have tiny impact on both the consensus 
time and the memory consumption of the agents. 

Examples of populations of mobile and communicat- 
ing agents pop out in many natural contexts, chiefly in 
cases of group of animals or micro-organisms. Yet, it 



is perhaps the technological advancement in the fields 
of robotics and telecommunications that makes the in- 
vestigation of this issue increasingly urgent. From this 
point of view, the results presented in this paper may be 
far-reaching. Both the mobility and the communication 
models we have adopted, indeed, are straightforward and 
might serve as a reference to gain insights into more com- 
plex and realistic models. At the same time, the scheme 
we have introduced can be extended in the future so as 
to address such issues as the coupling between motion 
and communication. For example, our results show that 
the larger the agents' mobility the shorter the consensus 
time, so that changing neighbors frequently turns out to 
be an efficient strategy. Thus, agents with a tendency to 
move apart from each other after a success could expedite 
convergence, preventing the formation of local clusters. 
On the contrary, individuals with a tendency to reinforce 
existing links by aligning their direction of motion after 
a successful interaction would probably lead to the emer- 
gence of different swarms or fiocks, internally agreeing on 
different conventions, that would move apart from each 
other and hinder the onset of a global agreement. Along 
the same lines, finally, another possibility for future work 
could imply the adoption and study of higher level com- 
munication protocols, such as spatially oriented language 
games [1, 
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